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ABSTRACT 

The problem of acquiring a simple but sufficiently accurate model of a dynamic system is made more difficult when 
the dynamic system of interest is a multibody system comprised of several components. A low order system model may 
be created by reducing the order of the component models and makmg use of various available multibody dynamics 
programs to assemble them into a system model. The difficulty is in choosing the reduced wder component modekto^ 
meet system level requirements. The projection and assembly method, proposed originally by Eke, solves this difficulty 
by forming the full order system model, performing model reduction at the system level using system level requirements, 
and then projecting the desired modes onto the components for component level model reduction. In this paper, the 
projection and assembly method is analyzed to show the conditions under which the desired modes are captured exactly-to 
the numerical precision of the algorithm. 

INTRODUCTION 

The problem to be solved is that of simulating the dynamics of a multibody system. A multibody system is 
comprised of two or more bodies or components connected at hinges. In general, the bodies may be rigid or fiexi e, 
and Ae hinges may have from one to six independent degrees of freedom. Often all deformations of each body from^ 
reference condition are in the linear range, although the resulting system dynamics is nonlinear. In this case nonlinear 
system models may be constructed using linear dynamic models for each component, but allowing large angle motion 
between components. This is the approach used in a number of existing multibody software tools. 

The problem is that system models constructed in this manner may be too large for use in control system design 
and simulation trades. Model reduction is needed to bring the model down to manageable size. If the system mode is 
available in linear form, system model reduction can be applied directly. For the class of inulubody problems 
discussed above, only the component models are available in linear form, and existing multibody software can be u 
if we reduce the component models before assembly into the system model. A multibody system is inherently a 
geometrically nonlinear system because of the time- varying, large-angle articulation between bodies. 

Component model reduction is typically done to some level anyway if the source of the model is a finite element 
program. This first level of model reduction often uses some simple criterion such as keep all cantilever modes below 
40 Hz." The challenge is to reduce the component model further in some manner that preserves how the component 
behaves when connected to the complete system; how the component affects system level requirements. The 
projection and assembly method described in this paper attempts to do this. 

Model reduction for linear systems has been addressed by a number of researchers, resulting in a variety of 
suggested linear system model reduction methods^. Fig . i gives a high-level view ofhow tltese . methods ^work 
Less attention has been paid to the problem of model reduction for components of multibody systems. Component 
modal synthesis methods 5678 have the capability of producing reduced order component models, but typically do so 
based on component-level rather than system-level criteria. When only one body in a mulubody system is flexib e, 
Macala 9 captures desired system modes exactly by augmenting the flexible body by the mass and inertia of the rigid 
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body. A subset of the free-free modes of this augmented body are then used as the flexible body component modes. 

Eke and Man 10 extend this capability to systems of more than one flexible body with a method that involves choosing 
system modes of interest, projecting the mode shapes of these desired modes onto each flexible component, reducing 
the order of each component accordingly, and assembling the components into a system model. Upon assembly, each 
of the original desired system modes is recovered exactly (to the numerical precision of the algorithm.) As can be seen 
in Fig. 2, this approach is conceptually more complicated than that shown in Fig 1., but allows the introduction of 
system level requirements. 
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Fig. 1. Conventional Model Reduction 



Fig. 2. Projection and Assembly Method 

This paper analyzes the method outlined in Ref. 10 to show why the desired modes are returned exactly, presents 
necessary conditions for the success of the procedure, and proposes an extension to the method to handle situations 
when these necessary conditions are not met. Simple examples are presented to demonstrate the workings of the 
algorithm. The name "Projection and Assembly Method" is used to describe this component model reduction method. 

DESCRIPTION OF METHOD 

The idea of the projection and assembly method is to decide what system modes are important and to choose 
component models which, when assembled, capture those important system modes. The projection and assembly 
method is described in detail in Ref. 10. It works as follows: 

•Acquire component models 

•Synthesize the system model in some configuration of interest 

•Apply any system level model reduction desired to choose which system free-free modes to retain 
•Project the mode shapes of these retained modes onto each component. 

•Choose new component states such that only these projected modes are admissible motions 
•Transform the component models into reduced order component models using these new states 
•Assemble the reduced order component models into a reduced order system model 
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The concept is made clearer by considering a simple qualitative example. Consider the planar motion of a system 
consisting of a rigid hub with two identical beam appendages, one on each side. Ignoring motions along the beam 
axes, the first five modes of this system are sketched in Fig. 3. 


Mode 1 





Fig. 3. Qualitative Beam Example 

In this example, the lowest three system modes are chosen to be retained in the reduced order model. Hie 
projection step is illustrated in Fig. 4; these system modes are projected onto each of the two components. The 
resulting projections are used as generalized mode shapes for component models. 
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Fig. 4. Example of Projection onto Components 

One side effect of the method is apparent by doing a little arithmetic. The projection and assembly method will project 
three modes onto each body. Assembling the components into a system gives two constraint relations (to match the 
halves of the rigid body together in angle and offset). When the reduced order component models are assembled, the 
reduced order system model will have four modes (3+3-2). These four include the three desired modes plus one 
"extraneous mode." 
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ANALYSIS 

Component Equations of Motion 


Assume we have /7 b bodies or components. The unconstrained equations of motion of each may be expressed as: 

Mj Xj+ KjXj= G'U, i = 1 , * • ■ , r) b 

where 

/ is the body index, 

Xj is a set of generalized coordinates describing the motion of body / as a free body in inertial 

space. This set of coordinates can be anything from geometric coordinates to free-free 
normal modes to cantilever modes augmented by six rigid body modes for the fixed 
end, 

Mj is the generalized mass matrix for body /, 

Kj is the generalized stiffness matrix for body /, 

U is the set of control inputs, 

Gj is the control distribution matrix for body /, and 
n b is the number of bodies. 

System Equations of Moti on , 

A multibody system is created by constraining the components to share certain common motions and by adding 
flexible connections between bodies. Assume that the constraints can be described in the form: 


AX = 0 2. 

where 

A=[ A. A 2 • * • A n ], X r = [ xjxl ■ • • x n J 

b 

Let n be the number of constraint equations in Eq. 2. The constraints may be introduced into the equations of 
c 

motion using a vector, A , of Lagrange multipliers. The constrained system is: 


Mj Xj + Kj Xj -GjU + aJ A, 

/= 1, • • • , n b 

3 

AX= 0. 


4 


Let P be any full rank matrix mapping a minimal system state, X, into X: 

X=Px\ or Xj= Pf(, i = 1 , • • ■ , rifr 5. 


The constraint equation becomes: 


APx= 0. 6. 

Since the states X are independent, AP= 0. Once P is chosen so that Eq. 6 is satisfied, the constraint equation (Eq. 4) 
is automatically satisfied. Inserting Eq. 5 into Eqs. 3 and pre-multiplying by P ; - gives: 
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/= V-'.'V 


7. 


T.T, 


Pj M'Pj X + p; K^PyX = P, G,u + Pj Aj A 


Summing over / : 


Where 


Mx+ Kx= Gu 

n b 

M= J^Uf, 

hi 



8 . 

9. 


G=XPfo 1 

hi 

Equation 8 is the system equation of motion incorporating all constraints. Converting Eq. 8 to modal form: 

x = &q 

q + G?q = <P t Gu . 


System Model Reduction 

Assume we choose some model reduction method which yields as its output a set of H R modes, q R to be 
retained with the remaining set of modes, q ^ , to be zeroed. Then we can partition 0 . 

*- [q”]- 

Or, setting 0. the reduced order system model is: 

5fl + ^ffi u 


jf = diag (coj) then a homogeneous solution to Eq. 12, and therefore also a solution to the system of equations 
3 & 4, is q = 9j COS( COjt). Each of Xj , Xp and A will similarly be described by sinusoids: 

Xj = Pj 0 R e j cos(cOjt), Xj = Pj& R ej {-o)f)cos{cofi, A = A Qj cos{(of). 

Inserting the above into Eq. 3 for L/=0 gives a relation which will be needed in a later derivation. 

+ K ; ]P, ^ = A\/ 16 ' 
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Component Model R eduction 

None of the above is unique to the projection and assembly method, which uses the above as a starttyg point. The 
concept is as follows: Cause each component to have, as an allowable motion, the mode shape of each retained mode 
projected onto the component. When the system is reassembled from reduced order components, the retained mode will 
still be an admissible lotion of the reduced order system. In the following, it will be shown that m addiuon to being 
an admissible motion of the reduced order system, it is a mode of the reduced order system. 


Consider the projection of Cf^ onto component /. Using Eqs. 5 and 15. 


*r PprflR 


17. 


In general, Cf^ should be of lower order than Xj. Where before, component / had H - degrees of freedom, Eq. 15 
restricts the motion to n R degrees of freedom. Let X Rj be a set of component / modes that span the space of 
component motions allowed by Eq. 15. In Ref. 10, the choice: X R j = q R is made, so: 

PptfH- 1£ 


Implicit in this choice is the assumption that the matrix Pf& R is of full column rank. This assumption is 
violated in a number of situations. The most obvious case is when one component has fewer degrees of freedom than 
the number of modes in <t> R Other examples arise when the projections of the modes are linearly dependant within 
the subspace of a particular component. In a later section of this paper, an alternative choice for X Rj is explored for 

situations where P® a is not of full rank. Writing the component equations of motion (Eq. 3) and constraint relation 
I hi 

(Eq. 4) in terms of the X Ri : 

M Ri x Ri + K Ri*Ri = G Ri U + A R ft’ ' = 1 > * n b 1 9 ' 

"b . 

A R*R ~ ^, A Ri X Ri = 

/= 1 

where ^ 

M„ K Ri K F p R 21 . 

G R<= 0 R P I a i A Ri- A PpR 22 ' 


and A n and Xr, are defined in the same manner as A and X . This system of equations in X Rj and A may be 
Hn »-> 


n t i • O 

formulated in terms of a minimal set of states, X ^ with some mapping r R . 


X R = P R X R' 

23. 

With this choice, Eq. 20 becomes: 

a r p r x r = 

24. 

Since the X R are independent: 
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A R P R - °- 


25 . 


In actual practice, P R has a specific form, but to understand the behavior of the reduced order system, we can consider 
any P^ which is of full rank and satisfies Eq. 25. If j = = 

desired retained modes, then Eq. 20 becomes: 


* * = X Dn , as will be the case for the 
Mn b 


( 


I Aft 


rd x Ri "° 


which is automatically satisfied in view of Eq. 6. This suggests that a partial choice for P ^ is the column: 

[ / / ■ • • / ] * . One full rank P R which satisfies Eq. 20 may be created by taking the singular value 
decomposition of a portion of A R . 


['W'W °i 


f A 1 


'A2 


-U# A1 V A1 


and choosing: 


p T - 

r R~ 


I 1 
0 


f A2 


so 


A R P R ~ 


[ a R 2 ’ • ’ A Rn ]v A2 


n b \ 

\ hi J 

' I'V’iVfl U A £ A 1 V T AI V A! 

V hi J 


as desired. Furthermore, P R is of full rank by construction. 

Starting from Eq. 19, the equations of motion in terms of X R are: 


= 0 


n b n b n b 

I pIMrPr:** + I p X/Vfl= I P^*pW^ 

hi hi hi 

The form of P R suggests a partitioning of X R and P R into desired and extra states: 


L X E. 


P nc ] 


P R st [ P RD ' RE\ 

ro 1 

where P RDi = I and P R £ “ V' J- 1° partitioned form, Eq. 28 is: 


26 . 


27 . 


28 . 
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"b 

i M Ri 

hi 


Z M Ri P REi 1- 


hl 


% T" ^ 


lL pT REi M Ri ^ \i P ReMrP RE i 
i=1 M 


EJ 


n*. 


Z^Ri Z'V REi 

j=1 i=1 


n b n b 

Yu pT REi K Ri Yj P ReFrP RE i 

hi 


b 

71 

mt 

hi 


' X D 1 

, x e\ 


"b 

hi 




iL pT REf*Ri 

hi 


u. 


29 . 


By construction, the system is capable of taking the shape of any of the n R desired modes. It remains to be shown 
that the Xq are free-free normal modes of the reduced order system. To show that they are requires only that 

X Q^jc 0 S 6)jt be a solution of Eq. 29 with U = 0. Assume X R = ^Q^jcOSCdyf, then: X R = ( ~0)j )X R and 
Eq. 29 becomes two equations: 

~ n b P n b 

Z w Ri ( _a> y) + Z^ffi 

hi hi 

'L P REi M R,(~ 0) f) + lL P REi K Ri 
hi hi 

If both left-hand sides in the above equations evaluate to zero, then the desired modes are modes of the reduced order 
system. Consider EM/?/ and 

n b n b n b 

j=1 '‘‘■'•hi hi • 

Similarly, 

Z^fl/ = ,fi2 ’ 

hi 

Eq. 30 becomes: 

((-a)jV + I2 2 )e / = = 0 

and so is satisfied. Consider Eq. 31: 


eposcojt = 0. 


31 


? _ 

eposcojt = 0 


30. 
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n b n b 

'L P JlEi M Ri(~ CO j) + X , P REpRi 

/= 1 


b 

71 

md 

hi 


e, = 

J (21) 


n u 




e i 


( = } X pT RE$fl A l A oj ( i) 'L P REi A Ri A 0 j - 1° V A2J*S\* (ii) 0A oj - 0 

hi 


Therefore the desired mode shapes and frequencies satisfy Eq. 29 and thus are normal modes of the reassembled reduced 
order system. 

Component Model Reduction — Extended Method 

As mentioned above, the choice: X Rj = q R depends on the matrix Pf& R being of full column rank. When this 
is not the case, the method can be extended to allow model reduction to proceed. Consider the singular value 
decomposition of suppressing the index i on the products in the SVD, let fbe the rank of Pf& R < l et Hr be 

the rank of O r , and let /7 ; . be the rank of P f (and the number of states in xp. If n=n^>n- (Ref. 10) 


PfP R .UIV T .[U t u 2 ] [ o']v r = u.iy. 

If r=np>n R (body / has few DOF) 


Pp R = UZV T = U[z 1 0] 


V 


= uz,vl 


If r<n f r<n R (linear dependant projected modes) 

Pp R = uzv T =[u 1 u 2 ] 


Z 1 0 

.0 0J 


= u,rv t , 


To ensure that the set Xpj is an independent set spanning the space of component motions, choose X Rj - 
Z^i) V^{i)qpf In the event that /(/)=n R , v](i) becomes V{i). This choice of X Rj gives for Eq. 18: 

x i~ U 1^ x Ri 

In the event that t{i)=n t Uj (/) becomes U{l). Define Q, = U^l). Eqs. 21-22 now take the form: 


M R rQ;Mp i 

G R r°! G i 


Km’°i K Pi 


A Ri~ A Pt 


Choosing: 


P T n- 


V,(1)Z,(1) V,(2)Z,(2) • • • v, (n b )z, (n b ) 


0 


'A 2 


18A. 

21A. 

22A. 

27A. 
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gives as desired, A^P^ = 0. Moreover, is again full rank by construction. Partition as before: 
P R = [ P RD P Re]' w^e P RDj = and P fle = [\£j. Eq. 29 becomes: 


n b n b 

'5L P RDi M Ri P RDi '£ P RDi M Ri P REi L 
i= i 


b 

h= 1 


n L 


n L 


^ P REi M Ri P RDi REi 

i=1 i=1 


l X E. 


n b n b 

Ys P RdPrP, RD i Ys P RdPrP RE i 
b=1 i=1 

n b n b T 

'£ P REi K Ri P RDi ^h P RE\^Ri P REi 

hi hi 


~ x d]_ 

. x e1~ 


% P RDi G Ri 

h 1 


n 


b 

► 

mJ 

/= 1 


'L P REi G Ri 


U 


29A. 


and the proof that the desired modes are normal modes of the reassembled reduced order system proceeds exactly as 
before, using the above definition of PpjQj and Eqs. 21 A and 22A. 


SIMPLE EXAMPLES 

One Dimensional Three Disk Example 



Fig. 5. One Dimensional Three Disk Example 


Consider Fig. 5. In this example, there are three disks, with rotational displacements (from left to right) y y y ^ 

and yj and inertias 4J, J, and J connected to ground and each other by torsion rods of equal spring constant, k. We 

choose to consider this simple system as being composed of two simpler subsystems of components. We divide the 
middle disk in half and allocate one half to each subsystem. Subsystem 1 contains the large disk and the left half of 

the middle disk. Take X [y y y^] ^ Subsystem 2 contains the rest of the system. Take X^= [y^ y$ ] Choosing 

units to make J and k equal to unity, the mass and stiffness matrices for each component are: 


Mi = 


4 

0 


0 

0.5 




0.5 0 

0 1 


K 2 = 


1 

-1 
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The constraint relation that connects the subsystems is that X ^(2) = X^l). Expressed in terms of a constraint matrix, 

A: 

A -l A, A z ], [ X, 4>|*J"° 

with 

=[o i ]. 4>=[- 1 ol 


One choice of P which reduces this to minimal form is: 

-1 0 0 

p_ 0 0.7071 0 

r ~ 0 0.7071 0 

0 0 1 


The system mass and stiffness matrices are: 



4 

0 

0 

, K= 

2 

0.7071 

o 1 

M = 

0 

0.5 

0 

0.7071 

1 

-0.7071 


L o 

0 

1 J 


L o 

-0.7071 

2 J 


The eigenvalue and eigenvector matrices for this system are: 


0.2803 0 0 

0 = 

0.4457 

-0.2153 

0.0703 

0 1.1694 0 

-0.5539 

-0.8155 

1.0140 

0 0 3.0502 J 


-0.2277 

-0.6943 

-0.6827 J 


System model reduction: Assume we wish to 


capture only the lowest frequency system mode (13^ = 0.2803), then 


Or = 


0.4457 

-0.5539 

L-0.2277J 


Component model reduction: choose X R1~ X R2~ CI R’ so 


X-i = Pi ®RX m = 


-0.4457 
L -0.3917 


*R1 x 2 = P2®R X R2 ~ 


-0.3917 

-0.2277 


X R2 


and the reduced order component mass and stiffness matrices are: 

MFn =[0.8714], Kr i =[0.2016], M m =[0.1286], /fe = [0.0787] 

The reduced order constraint matrix is: 

Ar = [ -0.3916 0.3916 L 


Choose P^ : 


Pr = 


-0.7071 
-0.7071 J- 


This gives the reduced order system: 

[0.5]Xfl + [0.1402]x fl = 0. 

Which has a single eigenvalue at = 0.2803. In this case, no extra modes are created because it happens that 
(2ripj- # constraints) = This is not true in general. The next example produces extra modes. 
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Fig. 6 . One Dimensional Five Disk Example 



Consider Fig. 6 . In this example, there are five disks, with displacements (from left to right) y ^y? Y ' 3 . 
y 4 , and y 5 and inertias 4J, J, J, J, and J connected to ground and each other by torsion rods of equal spring 

constant, k. We choose to consider this simple system as being composed of two simpler subsystems of components. 
We divide the middle disk in half and allocate one half to each subsystem. Subsystem 1 contains the large disk 
through the left half of the middle disk. Take X f = [/,, y z y 3 ] T . Subsystem 2 contains the rest of the system. Take 

Xg= [y\j, y^, y^] Choosing units to make J and k equal to unity, the mass and stiffness matrices for each 
component are: 


' 4 0 0 1 0.5 0 0 1 [ 2-101 [ 1 - 10 ' 

M\= 0 1 0 > m 2= 0 1 0 K i = -12-1. ^2= -12 -1 

. 0 0 0.5 J 0 0 1 J : Lo-llJ L0-12-- 

The constraint relation that connects the subsystems is that X ^(3) = 1)- Expressed in terms of a constraint 

matrix. A: 

A-lA, *>], [ = 

with 

=[ 0 0 1 ]. A >=[- 1 0 0]. 


One choice of P which reduces this to minimal form is: 

"-1 0 0 0 0 * 

0 -0.7071 0.7071 0 0 

p 0 0.5 0.5 0 0 

r = 0 0.5 0.5 0 0 

0 0 0 1 0 

0 0 0 0 1 J • 

The system mass and stiffness matrices are: 

"4 0 0 0 0 

0 0.75 -0.25 0 0 

M= 0 -0.25 0.75 0 0 

0 0 0 1 0 

L0 0 0 0 1 

The eigenvalue and eigenvector matrices for this system are: 


r 2 -0.7071 

-0.7071 2.2071 
K= 0.7071 -0.5 

0 -0.5 

L 0 0 


0.71 0 0 

-0.5 -0.5 0 

0.7929 -0.5 0 

-0.5 2 -1 

0 -1 2 J • 
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£2 


2 


0.1933 0 0 0 0 1 

0 0.5466 0 0 0 

0 0 1.4696 0 0 

0 0 0 2.6609 0 

. 0 0 0 0 3.6296J * 


0 = 


' -0.3489 0.3208 -0.1392 
0.1217 0.3655 -0.0438 
0.7270 0.4501 -0.8075 
0.3386 0.5329 0.3142 
L 0.1874 0.3666 0.5924 


-0.0704 

0.7626 

-0.0986 

0.3895 

-0.5894 


-0.0310 
0.8765 
0.3272 
-0.5924 
0.3635 J 


System model reduction: Assume we wish to capture only the two lowest frequency system modes {£2? = 0.1933, 
0.5466), then 

"-0.3489 0.3208 
0.1217 0.3655 
Or =\ 0.7270 0.4501 
0.3386 0.5329 
0.1874 0.3666-1- 


Component model reduction: choose x R1- X R2~ C IR' so 


X-, - P<\0RX m = 


0.3489 -0.3208 
0.4280 0.0598 
L0.4244 0.4078 J 


* 


*2 = P2&RXR2 = 


[0.4244 0.4078 
0.3386 0.5329 
LO.1874 0.3666J 


*R2 


and the reduced order component mass and stiffness matrices are: 


Mr i 


0.7602 -0.3357 

-0.3357 0.4985 


Kr] = 


0.1280 

-0.0831 


-0.0831 
0.3689 J’ 


Mr2 = 


0.2398 0.3357 
.0.3357 0.5015J’ 


K J 0.0653 0.0831 
0.0831 0.1777. ■ 


The reduced order constraint matrix is: 


Choose P R 


Ar = [0.4243 0.4078 -0.4243 -0.4078 ]- 


PR 


-0.8603 0 0 1 

0.2904 0.5926 0.5696 

-0.3021 0.7762 -0.2151 

-0.2904 -0.2151 0.7933 J- 


This gives the reduced order system: 


0.8954 0.1781 0.0875 ' 
0.1781 0.2307 0.2649 
. 0.0875 0.2649 0 . 3739 - 




" 0.2029 0.0883 0 . 0503 ' 
0.0883 0.1494 0.1383 
. 0.0503 0.1383 0 . 2062 - 


Xr = 


0 


which has three eigenvalues at £2? = (0.1933, 0.5466, 1.6873). The first two are the desired system modes, while the 
third does not match any of the original system modes; it is an "extraneous mode." 
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In this paper the model reduction method described by Eke and Man in Ref. 10 has been analyzed to 
demonstrate why the desired modes are returned exactly. An explicit set of necessary conditions involving the rank of 
the projection matrix has been presented, and an extension to the method has been proposed which removes those 
conditions. The method was demonstrated using two simple examples. 

Future work will address extending the method to handle variable configuration systems such as those with 
multiple articulation angles, better characterizing the "extraneous" modes which are a by-product of this method, and 
examining scaling issues which will arise when relative sizes of singular values are used to determine how many 
independent modes are projected onto a component 
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